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Abstract 

The bosonic IIB matrix model is studied using a numerical method. 
This model contains the bosonic part of the IIB matrix model conjectured 
to be a non-perturbative definition of the type IIB superstring theory. 
The large N scaling behavior of the model is shown performing a Monte 
Carlo simulation. The expectation value of the Wilson loop operator is 
measured and the string tension is estimated. The numerical results show 
the prescription of the double scaling limit. 
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1 Introduction 



Recent works in string theory have proposed some models as non-perturbative 
formulations Especially, the IIB matrix model 0,^] has been considered as 

a constructive definition of the type IIB superstring theory. This model is zero- 
volume limit of the ten-dimensional large iV supersymmetric Yang-Mills theory 
and defined by the following action, 

S = ~i* (\[A»,A v f + \i>T»[A^ , (1) 

where and ip are N x N traceless Hermitian matrices. The interesting feature 
of this model is that the space-time coordinates are considered as the eigenval- 
ues of these matrices. Then, we expect that the fundamental issues including 
the dimensionality and the quantum gravity can be understood by studying the 
dynamics of the model. 

To take the continuum limit (N — > oo) for the IIB Matrix model, a sensible 
double scaling limit should be determined dynamically. The scaling property 
of the model has been studied with the light-cone string field Hamiltonian of 
the type IIB superstring theory |4]]. The scaling property of the two important 
quantities of the model, the string scale (a') and the string coupling constant 
(g s tr) are determined as follows, 

a' ~ g 2 N a+b = g 2 N"< ~ Constant, 
g s tr ~ N a - b . (2) 

For the finite value of the string coupling constant (g s tr) > one have a restriction, 
a = b. In the IIB matrix model, the exponent (7) plays an important role that 
one can take the large N limit as the continuum limit for the IIB Matrix model. 
This exponent is determined dynamically. 

For studying the dynamical aspects of the model, some Monte Carlo simula- 
tions have been performed. In Ref. |5],|6|], the existence of the large N limit of 
bosonic SU(N) Yang-Mills matrix model for D > 2 has been discussed analyt- 
ically. Then, in Ref. 0, the bosonic model has been studied with a analytical 
method, the 1/D expansion, as well as a numerical one. The numerical results 
support the 1/D expansion as an effective tool detecting the large N scaling be- 
havior of the model. The leading term of the 1/D expansion at D > 3 suggests 
that the exponent (7) takes a value of 1. We study the model in ten dimensions 
and reconfirm the expected scaling behavior by a numerical method with a larger 
size matrix. 

In analogy with the two-dimensional model, the Eguchi-Kawai model, we 
study the scaling property of the Wilson loop in ten-dimensions. The area law 
of the Wilson loop operator has been found in the four-dimensional model in 
Ref. j§|. In this article, we also obtain that the area dependence of the Wilson 
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loop obeys the area law in ten dimensions. We calculate the string tension from 
the area law of the Wilson loop. We thus consider that the scaling property of 
the string tension is estimated in ten-dimensional model. 

This paper is organized as follows. In section 2, we review the model and some 
perturbative analysis. In section 3, we show the numerical results and the scaling 
property of the model. Then, we present the data which show the existence of 
the double scaling limit in the model. Finally, in section 4, we summarize and 
discuss our numerical results. 



2 Large N behavior of correlation functions 

First, let us remind the perturbative arguments of the bosonic model and describe 
briefly the large N behavior of the correlation functions of the gauge fields based 
on[0. 

The bosonic model of the IIB matrix model is given by 

S bosonic = ~ "~ gtrfAu, A v ] [A^, A u ], (3) 

where A n are N x N Hermitian matrices representing the ten-dimensional gauge 
fields. The coupling constant (g) is nothing but a scale parameter and is absorbed 
with the rescaling of the gauge field (A^) as A u — > -j=A n . 
The Schwinger-Dyson equation is given by 




leading to the relation of the correlation function, 

-<tT([A u ,A u } 2 )>=D(N 2 -l)g 2 . (5) 

For the estimation of the large N behavior of the correlation function, < 
tr([A M , A v ] 2 ) >, the matrices (A u ) are decomposed into the diagonal parts (X*) 

and the off diagonal parts (A u J ), and eq.(||) is taken up to the second order of 
the off diagonal elements (A^ 3 ), 

< tr[A M , A u ] 2 >~ 2< tr[X„, A V \[X^ A v \ > - < tr[X M , A U ][X U , A n ] >. (6) 

Counting the order of the diagonal parts, the large N behavior of the leading 
term of the correlation function counting with the order of the diagonal parts is 
given by 

<ti\A a ,A u ] 2 >~g 2 N 2 . (7) 

For the finite value of the correlation function, an upper limit of the typical scale 
of the extend of the space-time, for example can be R 2 =< trA 2 >, is suggested 
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as R 2 < N 1 ! 2 perturbatively and from the 1/D expansion the large N behavior 
is shown as |?J 

R 2 ~ gN 1/2 . (8) 

In the similar manner the correlation functions can be calculated perturbatively. 

Next, we consider the Wilson loop operator in the IIB matrix model. The 
Wilson loop operator and the large N behavior have been studied with the light- 
cone string field theory of the type IIB superstring The Wilson loop operator 
(w(C)) is defined as, 

w(C) = tr(v(C)), (9) 

where C denotes the closed path and v(C) is defined as v(C) = Up - ■ - U^ = 
Pa exp(z / dak^Afj) in the bosonic model. The matrices are considered as 
the unitary matrices, 

Up = exp(z / dlA^). (10) 



In the ordinary lattice gauge field theory, the expectation value of the Wilson 
loop operator which spreads a large area behaves as follows, 

w{I, J) ~ exp(—KI x J), (11) 

where I and J are the side lengths of the rectangular loop and K denotes the 
string tension. In the same analogy, we study the Wilson loop operator of the 
bosonic model of the IIB matrix model. 

From the scaling relation of two-dimensional Eguchi-Kawai model, 

g 2 N ~ Constant, (12) 

It is expected that the similar scaling relation also holds in the IIB matrix model 
as i 

a' ~ g 2 N^ ~ Constant. (13) 

The exponent (7) should be determined dynamically from the model. For the 
large N limit, the parameter {g 2 N" f ) must be fixed in the IIB matrix model in 
the same manners as the parameter (g 2 N) must be fixed in the Eguchi-Kawai 
model. 

We notice that the bosonic model is equivalent to the D > 2 Eguchi-Kawai 
model in the weak coupling limit. Since the U(1) D symmetry rotates all the eigen- 
values by the same angle, the following expansion is valid in the weak coupling 
region, 

Up ~ e^e a ", (14) 

where a M take constant values due to the U(1) D symmetry and are small. 
The bosonic model action can be obtained by expanding the action of the ten- 
dimensional Eguchi-Kawai model in terms of A^. When the higer order terms of 
A^ can be neglected, we can obtain the area law of the Wilson loop operator as 

w(I x J) =< tr(C^ • • • U^jxj >^< tr(e jA " • • • e^) IxJ >~ exp(-^(/ x J)). 

(15) 
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In Ref. ||, the area dependence of the Wilson loop operator has been measured 
and found the area law in D = 4. 

In following section, we will show that the area law holds in the bosonic model 
of the IIB matrix model in D = 10. 



3 Monte Carlo simulation and Large N behavior 
of the bosonic model 

For numerical simulation of the model, we consider the partition function, 

dAe~ Sbosonic , (16) 



where the action Sb osonic is given by eq.([|) and the measure of gauge fields is 
defined by 
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dA= n 



N N 

niW4?: 

i=i j=i 



(17) 



The action is quadratic with respect to each component, which means that 
we can update each component by generating gaussian random number in the 
heat-bath and the Metropolis algorithm. 

To confirm the large N behavior of the model, we measure the following 
expectation values, 

R 2 =< ^r(A 2 ) >, 
<±tr({A„A„] 2 )> . (18) 

From the perturbative calculation, the large N behaviors are shown as 

R 2 ~ gN l/2 , 

<±tr([A, 7 A v ] 2 )>~g 2 N. (19) 

In Fig.|l| and Fig||, we show the numerical results of the extent of space-time 
(R 2 ) and the correlation function (< tr([A M , A u ]) 2 >) for D — 10 with iV = 
16, 32, 48, 64, 128, respectively. We obtain the large N behavior as 

RZ-gN - 5 ^, 

<ltr([A M ,^] 2 )>~/iV- 1 - ( 1 ». (20) 

By the simulation using the larger size matrix, the numerical result get close to 
the expected results. 
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Figure 1: The measurement results of the extent of the space time. 



Then, we consider the Wilson loop operator, 



W (C) =< -^tr(e iA " 



■e iA ") > . 



(21) 



We take the loop (C) as the rectangular (J x J) where we select any two direction 
(/z, v) in ten dimensions. 

We show the measurement results of the loop operator in Fig.[3[ Since the 
dependence of the direction of the rectangular is not found, we consider that in 
the bosonic model the isotropy of the ten-dimensional space-time is not broken 
down spontaneously. 

Then, we calculate the large N behavior by the Wilson loop operator. From 
the numerical results, the Wilson loop operator, < jjtr(e lAfJ - ■ ■ ■ e tA ")j x j >, closes 
to the exponential curve with the large size area. It means that the Wilson loop 



operator obeys the area law eq.fllq). 

Then, we can obtain the string tension. In Fig.|4], we plot the string tension 

We find the large N behavior of the string tension (K), 

K~g- 2 N~ im ^. (22) 
From the string tension (K), the string scale (a') is estimated as 

a' ~ 1/K ~ 2 2 jV L07(1) = Constant. (23) 
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Figure 2: The measurement results of the extent of the correlation function of 

(40- 

We remark that the string field theory of the light-cone frame suggests |3[] 

a' ~ g 2 N = Constant. (24) 

We consider that the numerical result closes to the analytic one and that the 
four- dimensional model has the same scaling property ||. 

Furthermore the numerical result suggests that the large N behavior of the 
square root of the string tension approximately equal to the inverse of the extent 
of the space-time. 

K 1 / 2 ~ g^N" - 5 ^ „ R -\ (25) 

It means that the Planck scale of the theory has the same scaling property of the 
extent of the space time in ten-dimensions. It also holds on the two-dimensional 
model H and the four-dimensional model ||. 

4 Summary and Discussion 

Let us summarize the main points made in our calculation. We confirm that the 
Wilson loop operator in the ten-dimensional bosonic model obeys the area law 
similar to the two and four- dimensional model For the scaling behavior of 

the bosonic model of the IIB matrix model, our numerical estimation is 

a ~ g 2 N im{1) = Constant. (26) 
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Figure 3: The large area behavior of the Wilson loop operator w(I x J), 
dot corresponds to different choice of I and J. 
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Our results show the ten-dimensional space-time extends with g^jyo^io) 
and the extent scale of the space-time approximately equal to the Planck scale 
(I), I oj K~ x l 2 . Furthermore, these results are consistent to the suggestion from 
the string field theory on light-cone frame |3| and the 1/D expansion 0. From 
the numerical results, we consider that the ten-dimensional bosonic model and 
the four- dimensional model have the same scaling property. 

In this article, we calculate only the bosonic model of the IIB matrix model. 
For the future work, we are also considering the numerical simulation of the full 
model including the fermionic term. The supersymmetric four-dimensional model 
has been studied in Ref. ||. Optimistically, we expect that we can simplify the 
fermionic term with the perturbative calculation. In Ref. [EIJ, it is claimed that 
the model with the 1-loop effective action of the IIB matrix model produces the 
four- dimensional space-time from the ten-dimensional space-time. We thus make 
preparations the calculation of the modified model including the fermionic term. 
The improved supersymmetric model including the fermionic term is studied |[L1|| . 
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Figure 4: The measurement results of the string tension a. 
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